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Abstract A multi-purpose fixed-target experiment using the 
proton and lead-ion beams of the LHC was recently pro- 
posed by Brodsky, Fleuret, Hadjidakis and Lansberg, and 
here we concentrate our study on some issues related to the 
spin physics part of this project (refeiTed to as AFTER). We 
study the nucleon spin structure through pp and pd pro- 
cesses with a fixed-target experiment using the LHC pro- 
ton beams, for the kinematical region with 7 TeV proton 
beams at the energy in center-of-mass frame of two nucle- 
ons a/s = 115 GeV. We calculate and estimate the cos2(j) az- 
imuthal asymmetries of unpolarized pp and pd dilepton pro- 
duction processes in the Drell-Yan continuum region and at 
the Z-pole. We also calculate the sin(20 — (ps), sin(20 + (j)s) 
and sin 20 azimuthal asymmetries of pp and pd dilepton 
production processes with the target proton and deuteron 
longitudinally or transversally polarized in the Drell-Yan 
continuum region and around Z resonances region. We con- 
clude that it is feasible to measure these azimuthal asymme- 
tries, consequently the three-dimensional or transverse mo- 
mentum dependent parton distribution functions (3dPDFs or 
TMDs), at this new AFTER facility. 



1 Introduction 

Recently a multi-purpose fixed-target experiment using the 
proton and lead-ion beams of the Large Hadron Collider 
(LHC) extracted by a bent crystal, referred to as AFTER in 
the following, was proposed by Brodsky, Fleuret, Hadjidakis 
and Lansberg [1]. Such an extraction mode will not alter the 
performance of the collider experiments at the LHC. The 
center-of-mass energy is ^/snn =115 GeV with the LHC-7 
TeV proton beam and ^/snn = 72 GeV with the a lead run- 
ning with 2.76 TeV-per-nucleon beam, and it can be even 
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higher by using the Fermi motion of the nucleons in a nu- 
clear target. This project will provide a unique opportunity 
to study the nucleon partonic structure, spin physics, nuclear 
matter properties, deconfinement in heavy ion collisions, W 
and Z productions, exclusive, semi-exclusive and backward 
reactions, and even further potentialities of a high-energy 
fixed target set-up. We concentrate our study on some issues 
related to the spin physics for the AFTER proposal. 

The study of the three-dimensional or the intrinsic trans- 
verse momentum dependent distribution functions (3dPDFs 
or TMDs) has received much attention in recent years [2]. 
Such new quantities of the nucleon provide us a significant 
perspective on understanding the three-dimensional struc- 
ture of hadrons and the non-perturbative properties of quan- 
tum chromodynamics (QCD). The intrinsic transversal mo- 
mentum of partons may cause special effects in high energy 
scattering experiments [3]. Azimuthal asymmetries of unpo- 
larized and single polarized Drell-Yan processes are among 
the most challenging issues of QCD spin physics [4-6]. 

The first measurement of the angular distribution of Drell- 
Yan process, performed by NAIO Collaboration for 7zN, in- 
dicates a sizable cos 20 azimuthal asymmetry [7, 8] which 
cannot be described by leading and next-to-leading order 
perturbative QCD [9]. Furthermore, the violation of the Lam- 
Tung relation [10] which is obtained from the spin- 1/2 na- 
ture of quarks and the spin-1 nature of gluons, just like the 
Callan-Gross relation in the deep-inelastic scattering [11], 
was measured by Fermilab E615 Collaboration [12]. This 
violation was also tested by E866/NuSea Collaboration through 
the pd and pp Drell-Yan dimuon processes in recent year [13, 
14]. 

Large single spin asymmetries (SSAs) were observed 
experimentally in the process p/?^ -^ nX two decades ago [15- 
20]. SSAs in semi-inclusive deeply inelastic scattering (SIDIS) [21- 
27] with one colliding nucleon transversely polarized have 
also been measured by several experiments. Standard per- 



turbtive QCD based on collinear factorization to leading twist 
failed to explain these asymmetries [28]. 

The Drell-Yan process is an ideal ground for testing the 
perturbative QCD and for probing the 3dPDFs or TMDs, 
as it contains only the distribution functions with no frag- 
mentation functions, and its differential cross section is well 
described by next-to-leading order QCD calculations [29]. 
In this paper, we calculate azimuthal asymmetries of pp 
and pd dilepton production processes in Drell-Yan contin- 
uum region and around the Z-pole through a fixed-target ex- 
periment using the LHC proton beams with the proton or 
deuteron target unpolarized and transversally or longitudi- 
nally polarized. The paper is organized as follows. In Sect. 2 
and 3, we respectively calculate the azimuthal asymmetries 
in unpolarized and single polarized pp and pd processes. In 
Sect. 4, we present the numerical results of these asymme- 
tries. Then, a brief discussion and conclusion is contained in 
Sect. 5. 



2 The cos 20 azimuthal asymmetries of unpolarized pp 
and pd processes 

The Drell-Yan process is an ideal ground to investigate the 
hadron structure, because it only probes the parton distribu- 
tions without fragmentation functions. It was naively specu- 
lated that the polarization of at least one incoming hadron is 
necessary to investigate the spin-related structure and prop- 
erties of hadrons. However, it is not the case if we take the 
intrinsic transversal momentum of quarks inside the hadron 
into account. As mentioned before, the standard perturba- 
tive QCD to leading and next-to-leading order failed to de- 
scribe the sizable cos 20 azimuthal asymmetry and the Lam- 
Tung relation violation of the unpolarized Drell-Yan experi- 
ments [7, 8, 12-14]. Several attempts were made to interpret 
this asymmetry, such as the factorization breaking QCD vac- 
uum effect [9] (which corresponds possibly the helicity flip 
in the instanton model [30]), higher twist effect [31-33] and 
the coherent states [34]. Boer pointed out that the cos 20 az- 
imuthal asymmetry could be due to a non-vanished 3dPDF 
or TMD hf{x,pj) [35], named as the Boer-Mulders func- 
tion later, as one of the eight leading-twist 3dPDFs or TMD 
distribution functions contained in [36, 37] 
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The sizable cos 20 azimuthal asymmetry can arise from a 
product of two Boer-Mulders functions of two incoming 
hadrons by establishing a preferred transverse momentum 
direction from the spin-transverse momentum correlation. 
This effect is called the Boer-Mulders effect [35]. Many the- 
oretical and phenomenological studies are carried out along 
this direction [38-56]. 

The Boer-Mulders function hf, as well as the Sivers 
function f^, is a naively time-reversal odd (T-odd) distribu- 
tion function, characterizing the correlation between quark 
transversal momentum and quark transversal spin. There- 
fore, it was thought to be forbidden for a long time be- 
cause of the time-reversal invariance property of QCD [57]. 
However the model calculations taken by Brodsky, Hwang 
and Schmidt indicated that these non-vanished naively T- 
odd distribution functions, hf and f^ can arise from the 
final or initial state interaction between the struck quark and 
the target remnant in the SIDIS and Drell-Yan processes 
at leading-twist level [58, 59]. In general, the path-order 
Wilson Une, which arises from the requirement of a full 
QCD gauge invarient definition of 3dPDFs or TMD distribu- 
tion functions, provides non-trivial phases and leads to non- 
vanished T-odd distribution functions [60-65]. Due to the 
present of the Wilson line, opposite sign of Boer-Mulders 
function or Sivers function in SIDIS and Drell-Yan pro- 
cesses is expected [65-67] 
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This relation still awaits for experimental confirmation. For 
hadron productions in hadron-hadron collisions, the situa- 
tion is more involved, since colored objects exist in both 
the initial and the final states. The multiple final or initial 
state interactions will generate process-dependent 3dPDFs 
or TMDs which may be different from those in SIDIS or 
Drell-Yan process [68-71]. This is also viewed as the break- 
down of the generalized 3dPDF or TMD factorization in the 
inclusive hadro-production of hadrons [72]. 

The angular differential cross section for unpolarized Drell- 
Yan process has the general form: 
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where Q is the solid angle and A, fi, and v are angular dis- 
tribution coefficients. For azimuthal symmetrical scattering. 




Fig. 1 The CoUins-Soper frame. 

the coefficients ji ~v = 0. The polar and azimuthal angular 
and (j) are defined in the Collins-Soper (CS) frame [73], 
as shown in Fig. 1. It is the center of mass of the lepton 
pair with the z axis defined as the bisector of two incoming 
hadrons. The polar angular is defined as the angular of 
the positive lepton with respect to the z axis direction, and 
the azimuthal angular (j) is defined as the angular of the lep- 
ton plane with respect to the proton plane. In this frame the 
Lam-Tung relation is insensitive to the higher fixed-order 
perturbative QCD [74] or the QCD resummation [75-77]. 

Taking into account the Boer-Mulders distribution, we 
can express the unpolarized Drell-Yan cross section as 
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where h = qj/\qj\ is the direction of the transversal mo- 
mentum transfer, pj and kr are the transversal momentum 
of quarks in the nucleons, mm is the mass of a nucleon, a^,,, 
is the electromagnetic fine structure constant and e„ is the 
charge of the quark with the subscript a showing the flavor 
The structure function notation in the equation is defined as 
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(7) 
Then the cos 20 azimuthal asymmetry can be expressed as 
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If we take both y* and Z boson into account, the cross 
section is expressed as [35] 
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The coefficients Ki, K2, K^,, K4 are expressed as 
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where the combinations of the coupling constants are 
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and the Z boson propagator factors are expressed as 
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where 9w is the Weinberg angle. Then the cos 20 azimuthal 
asymmetry is 
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Another azimuthal dependent term is the sin 20 term in Eq. 
(9). However the sin20 term is XjQ^ suppressed. This sup- 
pression can be found from (13) and (17). 

For pd dilepton production processes, we assume the 
isospin relation. The distribution functions of m or m quark 
in proton is the same as those of t/ or J^ quark in neutron, 
and the distribution functions of d or d quark in proton is 
the same as those of m or S quark in neutron. We can also 
neglect the nuclear effect of deuteron, since it is a weakly 
bound state of a proton and a neutron. Therefore, for pd 
processes, we need to replace the 3dPDFs or TMDs of the 
target proton in Eq. (8)(18) as 
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with / representing f\ or h^ and similar for d, u and d 
quarks. Then we can get the cos 20 azimuthal asymmetry 
coefficient v for pd dilepton production in Drell-Yan con- 
tinuum region and around the Z pole. 



3 The sin(20 — 0s), sin(2(/) + ^s) and sin 20 azimuthal 
asymmetries of single polarized pp and pd processes 

Large SSAs observed experimentally [15-27] cannot be in- 
terpreted by the standard perturbative QCD based on collinear 
factorization to leading twist. As a challenging issue in hadron 
structure and QCD spin physics, many theoretical studies 
were proposed to explain origin of such asymmetries [78- 
83]. In the 3dPDF or TMD framework, the non-vanished 
naively T-odd Sivers function f^ in Eq.(l), which char- 
acterizes the correlation between quark transversal momen- 
tum and hadron transversal spin, was applied to explain the 
SSAs observed in the process pp^ —>■ nX [79-82]. SSAs 
contributed by this Sivers effect in SIDIS processes with 
one nucleon transversally polarized have been measured by 
several experiments in recent years [21-23, 26, 27, 84-86]. 
The data on the Sivers SSAs have been utilized by differ- 
ent groups to extract the Sivers function of the proton on the 
basis of the 3dPDF or TMD factorization [67, 87-93]. 

For a fixed-target experiment, it is convenient to polar- 
ized the target to allow the SSAs measurements. Five lead- 
ing twist 3dPDFs or TMDs in Eq. (1), /j^, hf, h{j, gir and 
/ij^, vanish upon integrating over the transversal momen- 
tum kj- The two naively T-odd distribution function f^, the 
Sivers function, and hf, the Boer-Mulders function are ac- 
count for the SSAs in various processes. Four leading twist 
3dPDFs or TMDs, /iij-, /i|, hjj and /i|^, are chirally odd, 
so they describe densities of the probed quarks with helicity 
flipped. The hir and lijj have the relation with the hi that 
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Fig. 2 The definition of azimuthal angles. 



Q^ region. Then the leading order of the differential cross 
section can be expressed as [35, 108] 
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The distribution functions hi and hf respectively charac- 
terize the densities of transversely polarized quarks inside 
a transversely polarized proton and an unpolarized proton. 
The distribution functions hfj and hf^, arising from the dou- 
ble spin correlations in the parton distribution functions, re- 
spectively describe the densities of transversely polarized 
quarks in a transversely orthogonally polarized proton and 
longitudinally polarized proton. 

The chiral-odd 3dPDFs or TMDs are rather difficult to 
be probed in high energy scattering experiments, because 
they only manifest their effects by combining with another 
chiral-odd function, Collins fragmentation function in SIDIS 
or another chiral-odd distribution function in Drell-Yan. Some 
efforts have been made to extract the transversity from SIDIS 
data [94, 95] and to extract the Boer-Mulders function from 
SIDIS and Drell-Yan data [47, 49, 50, 96]. There are some 
extensive model calculations of hfj and hf^ [83, 97-107]. 

If the transverse momentum of the dilepton in the Drell- 
Yan process qj, is measured, we can apply the 3dPDF or 
TMD factorization [62, 67, 92, 93], which is valid in q^ < 



where the structure function notation is defined as Eq. (7), 
and the azimuthal angles and 0s are defined as shown in 
Fig. 2. Therefore, one can define the following azimuthal 
asymmetries: 
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If we take both 7* and Z boson into account, the cross 
section is expressed as [35] 
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where the coefficients Ki{9), ^3(6) and /r4(0) are defined 
in (10)-(13). Then the azimuthal asymmetries defined in 
(22)-(24) are expressed with Z taken into account as 
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The cos 20, cos (20 — 0s) and cos (20 + 0s) azimuthal de- 
pendent terms are 1/Q^ suppressed. This suppression can 
be found from (13) and (17). 

For pd dilepton production processes with the deuteron 
longitudinally or transversely polarized, we assume the isospin 
symmetry and neglect the nuclear effect as we do for unpo- 
larized pd processes. Therefore, we can get the sin (20 — 
0s), sin (20 + 0s) and sin 20 azimuthal asymmetries of pd 
dilepton production processes in Drell-Yan continuum re- 
gion and around the Z pole by replacing the distribution 
functions of the target via (19). 

4 Numerical results 

In this section, we calculate the azimuthal asymmetries of 
the pp and pd dilepton production processes with the proton 
or deuteron target unpolarized and longitudinally or trans ver- 
sally polarized in Drell-Yan continuum region and Z reso- 
nance region respectively. We present a numerical estima- 
tion of these azimuthal asymmetries for measurement at a 



fixed-target experiment using the LHC beams proposed by 
Brodsky, Fleuret, Hadjidakis and Lansberg [1]. With the 7 
TeV proton beams, the center-of-mass frame energy ^/s = 
115 GeV for two nucleons. 

The cross section of Drell-Yan process can also be ex- 
pressed depending on y and Q- instead of jci and X2 with just 
a Jacobian multiplied as 
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Xl is the momentum fraction of the parton in the beam pro- 
ton, and X2 is that of the parton in the target nucleon. For the 
single polarized processes, X2 is the momentum fraction of 
the parton in the polarized target nucleon, sometimes labeled 
as x^ in the literature. To calculate the azimuthal asymme- 
tries depending on Q, xp or qj, we need to integrate over 
the other variables in the numerator and the dominator of 
the expression of the asymmetries respectively. The rapid- 
ity is cut in [—4.8, 1] which is the easiest region to carry on 
measurements as discussed in [1], and this is where the mo- 
mentum fraction of the parton inside the polarized nucleons 
is the largest. 

In our calculation, we adopt the Boer-Mulders function 
h\ extracted from the unpolarized pd and pp Drell-Yan 
data [13, 14, 47, 50]. The parametrization of h^ for both 
valence and sea quarks has the form: 
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where the subscript "bm" stands for the Boer-Mulders func- 
tions, and q = m, d, u and d. The possible range of the pa- 
rameters Hg allowed by the positivity bound [109] can be 
described by a coefficient o which balance the contribu- 
tions of quark and antiquark. Hq — > (oHq for q = u, d and 
Hq -^ co^^Hq for q ^ u, d will not change the calculated 
cos 20 asymmetry in the unpolarized pp and pd Drell-Yan 
data. The range of (0 is 0.48 < (0 < 2.1, and we choose the 
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Fig. 3 The cos 2(j) azimuthal asymmetry depending on Q of unpolar- 
ized pp Drell-Yan process with both 7* and Z taken into account and 
allowed rapidity integrated in the cut [—4.8, 1]. The same cut of rapid- 
ity is chosen in Figs. 7, 11-16. 
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Fig. 4 The cos 2(j) azimuthal asymmetry depending on Xf of unpolar- 
ized pp Drell-Yan process at Q = 2 GeV. 



case = 1, which corresponds to the central values of H^, 
in our calculation. 

In Figs. 3 and 7, we show the cos 2^ azimuthal asym- 
metry depending on Q from 2 GeV to 30 GeV of the unpo- 
larized pp and pd Drell-Yan process at AFTER including 
Z taken into account. Figs. 4, 5, 8 and 9 respectively show 
the cos 20 azimuthal asymmetry depending onxf of the un- 
polarized pp and pd processes with Q~2 GeV and Q = 5 
GeV as for low and mid Q Drell-Yan regions at AFTER. 
Figs. 6 and 10 respectively show this azimuthal asymmetry 
of pp and pd processes at the Z pole at AFTER. 

To calculate the SS As of the pp and pd dilepton produc- 
tion processes at AFTER, we also need the distribution func- 
tions hfj,, hf^ and hi besides the Boer-Mulders function hf. 
In our calculation, we adopt ansatz of these T-even distribu- 
tion functions calculated from the light-cone quark-diquark 
model. In this model, the Melosh-Wigner rotation, which is 
important to understand the proton spin puzzle due to the 



Fig. 5 The cos2(|) azimuthal asymmetry depending on Xf of unpolar- 
ized pp Drell-Yan process at g = 5 GeV. 
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Fig. 6 The cos 2(p azimuthal asymmetry depending on qr of unpolar- 
ized pp process in Z resonance region. 



relativistic effect of quark transversal motions [110-113], is 
taken into account. This model has been applied to calculate 
helicity distributions [114], transversity distributions [115, 
1 16] and some other 3dPDFs or TMDs [100, 104, 1 17], and 
has been used to analysis related azimuthal spin asymme- 
tries in SIDIS processes [118, 119]. The model results of 
these distribution functions are expressed as [ 100, 104, 1 15, 
116] 



-i/;'^(x,4)w^(x,4), (34) 

J2x,kl) = - ^-f,,{x,kl)wl{x,kl), (35) 



where j = hi, hjj, hf^, and the superscript v stands for va- 
lence quark distributions. The factors W^,y{x,kj) are the 
Melosh-Wigner rotations for scalar or axial vector spectator- 
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Fig. 7 The cos 2(j) azimuthal asymmetry depending on Q of unpolar- 
ized pd Drell-Yan process with both y* and Z talcen into account. 



Fig. 9 The cos2(^ azimuthal asymmetry depending on Xf of unpolar- 
ized pd Drell-Yan process at Q = 5 GeV. 
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Fig. 8 The cos 2(/) azimuthal asymmetry depending on Xf of unpolar- 
ized pd Drell-Yan process at 2 = 2 GeV. 



Fig. 10 The cos 2<p azimuthal asymmetry depending on qr of unpolar- 
ized pd process in Z resonance region. 



diquark respectively, having the forms: 



W^'(x,*2) 



W>(X,4): 



<'-(x,4): 
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The distribution functions h\, lifj, lifj^ of sea quarks are 
constrained by the positivity bounds [109]: 



hlg{x,kl) <flg{x,kl), 



2m 



^iiinM,4: 






<flq{^,kj), 



<figix,kl). 



(40) 



(41) 



(42) 



When considering the effects of these distribution functions 
of sea quarks contributing to the asymmetries, we can get 
the upper and lower limits of the azimuthal asymmetries by 
saturating the positivity bounds. 

In Figs. 11-13, we respectively show the sin(20 — ^5), 
sin(2^ + ^s) and sin20 azimuthal asymmetries depending 
on Q from 2 GeV to 30 GeV of the target proton polarized 
pp Drell-Yan process at AFTER including Z taken into ac- 
count, and the corresponding results for deuteron target po- 
larized pd process are shown in Figs. 14-16. These asymme- 
tries in the far backward region, with rapidity cut [—4.8,-2], 




Q(GeV) 

Fig. 11 The sm{2(j) — (ps) azimuthal asymmetry Ajy depending 

on Q of target proton polarized pp Drell-Yan process witli both y* 
and Z taken into account. The dashed curves show the range of the 
asymmetry by considering the additional distributions of sea quarks 
constrained by the positivity bounds, corresponding to the same case 
as Figs. 18^0. 




Q(GeV) 

Fig. 12 The sm(2(j) + (ps) azimuthal asymmetry A"'^ depending 

on Q of target proton polarized pp Drell-Yan process with both 7* and 
Z taken into account. 
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Q(GeV) 

Fig. 13 The sin2<^ azimuthal asymmetry A^™ depending on Q of 
target proton polarized pp Drell-Yan process with both 7* and Z taken 
into account. 
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Q(GeV) 

Fig. 14 The sm{2(j> — (j>s) azimuthal asymmetry A™^ depending 

on Q of target deuteron polarized pd Drell-Yan process with both y* 
and Z taken into account. 




Q(GeV) 

Fig. 15 The sm(2(j) + (ps) azimuthal asymmetry A^J, depending 

on Q of target deuteron polarized pd Drell-Yan process with both 7* 
and Z taken into account. 
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Fig. 16 The sm2(j) azimuthal asymmetry A'^^ depending on Q of 
target deuteron polarized pd Drell-Yan process with both 7* and Z 
taken into account. 
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Q(GeV) 

Fig. 17 The sm{2(j) — (ps) azimuthal asymmetry Aj'^ depending 

on Q of target proton polarized pp Drell-Yan process with both 7* 
and Z taken into account and allowed rapidity integrated in the cut 
[—4.8, —2]. The same cut of rapidity is chosen in Figs. 18-22. 




Q(GeV) 

Fig. 18 The sm{2(j> + (j>s) azimuthal asymmetry A^y depending 

on Q of target proton polarized pp Drell-Yan process with both 7* and 
Z taken into account. 




Fig. 19 The sm2(j> azimuthal asymmetry A^" depending on Q of 
target proton polarized pp Drell-Yan process with both 7* and Z taken 
into account. 




Fig. 20 The sm{2(j> — (j>s) azimuthal asymmetry A™^ depending 

on Q of target deuteron polarized pd Drell-Yan process with both y* 
and Z taken into account. 
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Fig. 21 The sm{2(j) + (ps) azimuthal asymmetry A^J, depending 

on Q of target deuteron polarized pd Drell-Yan process with both 7* 
and Z taken into account. 
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Fig. 22 The sm2(j) azimuthal asymmetry A'^^ depending on Q of 
target deuteron polarized pd Drell-Yan process with both 7* and Z 
taken into account. 
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Fig. 23 The sin(2(^ — ^s) azimuthal asymmetry A"^ depending Fig. 26 The sin(2(|) — ^s) azimuthal asymmetry A™^ depending 

on Xf of target proton polarized pp Drell-Yan process at Q = 2 GeV. on Xf of target deuteron polarized pd Drell-Yan process at Q = 2 GeV. 
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Fig. 24 The sin(2(|) + <|)s) azimuthal asymmetry Ajjj depending Fig. 27 The sin(20 -f 0s) azimuthal asymmetry Aj-^ depending 

on Xf of target proton polarized pp Drell-Yan process at Q = 2 GeV. on Xf of target deuteron polarized pd Drell-Yan process at (9 = 2 GeV. 
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Fig. 25 The sin 20 azimuthal asymmetry A^" depending on Xf of Fig. 28 The sin 20 azimuthal asymmetry A^" depending on Xf of 



target proton polarized pp Drell-Yan process at Q = 2 GeV. 



target deuteron polarized pd Drell-Yan process at Q = 2 GeV. 
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Fig. 29 The sm{2(j) — (ps) azimuthal asymmetry Aj'^ depending Fig. 32 The sm{2(j> — (ps) azimuthal asymmetry A™^ depending 

on Xf of target proton polarized pp Drell-Yan process at Q = 5 GeV. on Xf of target deuteron polarized pd Drell-Yan process at Q = 5 GeV. 




Fig. 30 The sm{2(j) + (ps) azimuthal asymmetry A 
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depending Fig. 33 The sin(20 + 0s) azimuthal asymmetry A' 



sm{2(()+fc) 



depending 



on Xf of target proton polarized pp Drell-Yan process at Q = 5 GeV. on Xf of target deuteron polarized pd Drell-Yan process at Q = 5 GeV. 
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Fig. 31 The sin 20 azimuthal asymmetry A^" depending on Xf of Fig. 34 The sin 20 azimuthal asymmetry A^" depending on Xf of 



target proton polarized pp Drell-Yan process at Q = 5 GeV. 



target deuteron polarized pd Drell-Yan process at 2 = 5 GeV. 
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Fig. 35 The sm{2(j) — (ps) azimuthal asymmetry Ajy depending Fig. 38 The sm{2(j> — (j>s) azimuthal asymmetry A™^ depending 

on qr of target proton polarized pp dilepton production process at the on qr of target deuteron polarized pd dilepton production process at 
Z pole. the Z pole. 
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Fig. 36 The sin(2(|) + (pg) azimuthal asymmetry A^^ depending Fig. 39 The sin(20 + 0s) azimuthal asymmetry A^J^ depending 

on qr of target proton polarized pp dilepton production process at the on qr of target deuteron polarized pp dilepton production process at 
Z pole. the Z pole. 
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Fig. 37 The sin 20 azimuthal asymmetry A^^ depending on qr of F'g- 40 The sin 20 azimuthal asymmetry A^JJ, depending on qr of 



target proton polarized pp dilepton production process at the Z pole. 



target deuteron polarized pd dilepton production process at the Z pole. 
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are shown in Figs. 17-22, with Q running from 1 GeV to 10 
GeV. In Figs. 23-34, we respectively show the sin(20 — ^5), 
sin (20 + 05 ) and sin 20 azimuthal asymmetries depending 
on Xf of the target proton and deuteron polarized pp and pd 
Drell-Yan processes at AFTER with Q = 2 GeV and Q = 5 
GeV as for low and mid Q regions. Figs. 35-37 show the 
sin (20 — 0s), sin (2 + 05) and sin 20 azimuthal asymme- 
tries of the target proton polarized pp process around the 
Z-pole at AFTER. The corresponding azimuthal asymme- 
tries of pd processes with target deuteron transversally or 
longitudinally polarized are shown in Figs. 38-40. 



5 Discussion and conclusions 

In this paper, we calculate the cos 20 azimuthal asymmetries 
of unpolarized pp and pd dilepton production processes in 
the Drell-Yan continuum region and around the Z resonance 
region. We also calculate the sin(20 — 0^), sin(20 + 0^) and 
sin 20 azimuthal asymmetries of single transversally or lon- 
gitudinally polarized pp and pd dilepton production pro- 
cesses in these regions. 

Our calculations are concentrated on some issues related 
to the spin physics part of the AFTER project, a multi-purpose 
fixed-target experiment using the proton and lead-ion beams 
of the LHC extracted by a bent crystal, proposed by Brod- 
sky, Fleuret, Hadjidakis and Lansberg [1]. We present an es- 
timation of the azimuthal asymmetries for a fixed-target ex- 
periment using the LHC 7 Te V proton beams with the proton 
or deuteron target unpolarized and transversally or longitu- 
dinally polarized. As the target is conveniently polarized, 
it is an ideal ground to study the spin physics at AFTER 
with v^ =115 GeV and high luminosity. It is feasible to 
measure these azimuthal asymmetries at AFTER. This will 
help us to study the three dimensional or transverse momen- 
tum dependent parton distributions (3dPDFs or TMDs), and 
consequently help understand and test the QCD and hadron 
structure at such a high laboratory energy. 
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